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Abstract 

The reduction of a continuous Markov process with muhiple metastable states to a discrete rate 
process is investigated in the presence of slow time dependent parameters such as periodic exter- 
nal forces or slowly fluctuating barrier heights. A quantitative criterion is provided under which 
condition a kinetic description with time dependent frozen rates applies. Finally it is shown how 
the long time behavior of the underlying continuous process can be retrieved from the knowledge 
of the discrete process by means of an appropriate random decoration of the discrete states. As 
a particular example of the presented theory an over-damped bistable Brownian oscillator with 
periodic driving is discussed. 
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I. INTRODUCTION 



The separation of time scales is a frequently met dynamical feature of physical, chemical 
and biological processes. In practically all cases the parameters characterizing a particular 
system are only constant on a limited time scale. On longer time scales they may change 
their values because the surrounding environment slowly changes its state. We just note 
that even the most fundamental "constants" of nature like e.g. the fine structure constant 
possibly undergo a very slow change in time [1]. On the other hand, most natural processes 
also display very fast changes which usually will be disregarded in a phenomenological de- 
scription. They only will show up in a microscopic description in which the effective degrees 
of freedom of the phenomenological model are represented as functions of microscopic de- 
grees of freedom. The time evolution of these microscopic degrees of freedom typically takes 
place on much faster time scales than those of the macroscopic observables whereas, as al- 
ready mentioned, global aspects of the environment may undergo temporal variations that 
are much slower than the dynamics of the considered macroscopic system. 

It is this separation of time scales that allows one to describe a system in terms of a 
comparatively simple model that is closed in the sense that only a few macroscopic observ- 
ables are subject to the dynamics. On the slow time scale the key technical notion here is 
that of a constraint equilibrium that is reached by the fast part of a system before the slow 
constituents change their values on the slower time scale. This constraint equilibrium then 
immediately adapts to the changes of the slow system. 

This approach has been applied in many branches of physics and other sciences and is 
vital for various methods like adiabatic elimination, averaging, synergetics, subdynamics, 
chemical kinetics and hierarchical computational schemes [2]. In the present paper we con- 
sider Markovian processes the dynamics of which can be subdivided into a fast and a slow 
component. In particular, we have in mind the case where the slow motion is caused by 
transitions between metastable states [3] . The relaxation toward these states constitutes the 
fast part of the dynamics. Additionally, some of the system's parameters may undergo slow 
variations in time that may be caused by a slow drift, or fluctuations of the parameters, or 
by an external driving of the system. The dynamics of the considered system then has two 
slow components: The transitions and the driving. About their relation we do not make 
further assumptions: The slow intrinsic dynamics can be faster or slower than the external 
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driving that must only be slow compared to the fast dynamics of the system. In Ref. [4] we 
have called this type of situation the semiadiabatic limit, in contrast to the adiabatic limit 
in which the driving is slow compared to the total intrinsic dynamics including that of the 
transitions. Clearly, the adiabatic limit is contained in the semiadiabatic limit as a special 
case. 

Stochastic resonance is just one example for which the semiadiabatic limit is of relevance 
when the period of the signal is large [5, 6]. Ratchets, or Brownian motors, provide another 
class of systems which also operate in the semiadiabatic regime [7] . Atomic force microscopes 
and optical tweezers often act on time scales that are slow compared to the vibrational and 
other intrinsic time-scales of the molecules that are manipulated by them [8]. 

The paper is organized as follows. In Sect. II the reduction of a Fokkcr-Planck process 
with multiple mestastable states to a master equation is reviewed as it was developed in 
Refs. [9, 10]. This method is generalized to processes with slowly time dependent parameters 
in Sect. III. The time rate of change of the probability now consists of two contributions. 
One is given by the rate as it would result if the parameters of the process were frozen and 
the other contribution mainly takes into account the changes of the geometry of the domains 
of attraction of the mctastable states. This geometric contribution is proportional to the 
time rate of change of the parameters and therefore is negligible compared to the frozen 
rates if the time dependent parameters change sufficiently slowly. In contrast to previous 
investigations, [11] we find that the fastest time dependence of the parameters for which one 
may neglect the geometric contributions is not only determined by the deterministic time 
scales of the process but in general also depends on the strength of the noise. The way how 
the noise here enters depends on the particular process. 

The present investigation is complementary to previous work that is restricted to periodic 
driving with an intermediate regime of external driving frequencies [12]. 

In Sect. IV we assume the validity of the master equation with frozen rates and determine 
a time dependent decoration of the discrete states such that the dynamics of the continuous 
variables of the underlying Fokker-Planck process are recovered on the long time scale. The 
archetypical example for stochastic resonance namely that of a periodically damped bistable 
Brownian oscillator is discussed in Sect. V and the validity of the McNamara-Wiescnfeld 
model of stochastic resonance [13] is discussed. The paper ends with a summary in Sect. VI. 
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II. MULTISTABLE SYSTEMS AT SMALL NOISE 



In this section we review the reduction of a Fokker-Planck equation describing a time- 
homogeneous process with multiple metastable states to a master equation that gives the 
dynamics on the time scales of the transitions between the metastable states. For this 
purpose we consider a dynamical system under the influence of weak random forces and 
assume that these forces can be modeled by Gaussian white noise. Then the stochastic 
dynamics is described by a Fokker-Planck equation [14]: 

^p(x,t)=Lp(x,i) (1) 

where L denotes the Fokker-Planck operator: 

i--E|-^.W + Ea5;-A,(x) (2) 

I l,J 

Here, x = {xi,X2, denotes a point with coordinates Xi in the n-dimensional state space 
E; further, K(x) = (A'i(x)) is the drift vector and D(x) = (Z}jj(x)) the diffusion matrix 
resulting from the random forces. We further restrict ourselves to systems with a uniquely 
defined stationary probability density po(x) with respect to which detailed balance holds. 
Consequently, the Fokker-Planck operator satisfies the relation [14]: 

Lpo = poL+. (3) 

where L"*" is the backward, or adjoint Fokker-Planck operator: 

The tilde denotes the operation of time reversal, i.e. Xi — eiXi with parity = ±1 and po 
is the multiplication operator with the probabihty density, i.e. Po/(x) = po(x)/(x) where 
/(x) is an arbitrary state space function. 

In the deterministic limit the diffusion matrix goes to zero and the drift vector approaches 
the deterministic vector field K^^'' (x) governing the deterministic motion of the system: 

x(t)=K(°)(x(t)) (5) 

In the presence of weak random perturbations K(x) may differ from K°(x) by small noise 
induced contributions. 
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FIG. 1: The overdamped deterministic motion x{t) = —V'{x{t)) in the above shown potential 
induces a partitioning of the state space S = Mi into the domains of attraction Da of the three 
attractors a = 1,2,3 at the potential minima. Here the locations of the local potential maxima 
coincide with the boundaries between the domains of attraction. 

We here are interested in cases where the deterministic system (5) has a number m of 
coexisting attractors labeled by a. To each attractor a there belongs a domain of attraction 
T>a- The domains are disjoint and partition the total available state space: 



For weak noise the stationary probability density is almost zero everywhere except close 
to the attractors. If the noise vanishes the stationary density shrinks to Dirac 5-functions 
sitting at the attractors of the deterministic system. Because no transitions between the 
different attractors may occur in the noiseless case, each of the delta functions is a sta- 
tionary solution of the Fokker-Planck equation d / dxiKf {'x) po{x.) = belonging to a m-fold 
degenerate eigenvalue 0. 

The degeneracy also follows from the backward equation which is governed by the adjoint 




(6) 



m 




(7) 



a=l 
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Fokker-Planck operator L+ defined in eq. (4). In the deterministic case the diffusion matrix 
vanishes, so the second derivatives disappear and the stationary solutions of the backward 
operator are constant along the trajectories of the deterministic system. Because within a 
domain of attraction all trajectories approach the same attractor the stationary solutions of 
the noiseless backward equation are constant on the domain of attraction T>a- One then can 
choose the characteristic functions Xa\^) of the domains of attraction as m independent 
solutions of the stationary deterministic backward equation Kl°\-x.)dx''^\'x.) / dxi — 0: 



In the presence of weak noise these functions can be modified slightly such that they still 
solve the stationary backward equation up to corrections that are exponentially small in 
the noise strength: The step-like discontinuities then are smeared out on layers located 
at the boundaries of the domains of attraction Pq,. The functions can be constructed by 
means of the method of matched asymptotic expansions [15], see also Refs. [9, 10]. Outside 
the boundary layers no modification is required and hence they retain their property of 
localizing the different domains of attraction. In the following, these localizing functions 
will be denoted by Xa(x)- We summarize their main properties: They solve L+Xa(x) = 
up to exponentially small corrections in the noise strength, have the values one within and 
zero outside and interpolate between these values smoothly on a thin boundary layer 
[16, 17]. As solutions of these particular boundary value problems, the functions ^^(x) are 
the splitting probabilities that give the relative frequencies with which a trajectory starting 
at X first reaches the attractor a before it visits any other attractor [18]. 

The splitting probabilities in general are not eigenfunctions of the backward operator. 
However, they span the subspace of eigenfunctions that emerges from the eigenspace of the 
m-fold degenerate eigenvalue of the deterministic dynamics if weak noise perturbes the 
system and lifts the degeneracy. The resulting finite eigenvalues are still small compared 
to all other finite eigenvalues. Hence, this subspace describes the slow dynamics of the 
system at weak noise. In particular, it contains the constant function which is the exact 
eigenf unction of the backward equation belonging to the eigenvalue 0. At weak noise this 
eigenfunction can be represented as the sum over the characteristic functions Xa(x): 




for X G P, 



else 



a 



(8) 
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Because of detailed balance, the slow subspace of the forward operator is spanned by 
the basis set {xa(x)po(x)}, where Xa(x) = Xa(x) denotes the image of Xa(x) under time 
reversal. Hence, the probability density describing the slow dynamics can be represented as: 



m 



t) = 2_^ c«(i)xa(x)po(x) (10) 

a=l 

As already noted the characteristic functions Xa;(x) and, therefore also their time reversed 
partners Xq:(x), are well localized functions as long as the noise does not become too large. 

As a consequence, one may use the functions Xa(^) order to determine the probability 
Pait) that the system resides in the domain of attraction of the attractor a [10]: 

Pa{t) = {Xa.p{t)) (11) 

where the scalar product is defined as the integral of the product of its arguments over the 

state space: 

(/,^) = ^d^x/(x)^(x). (12) 
If we assume that the probability density p{t) is of the form (10) we obtain 

Pa{t) = Ca{t)na (13) 

where we have made use of the localization of the characteristic functions at weak noise: 

(Xa,Xa'Po) = 5a,a'na (14) 

and have introduced the function Ua = J dxXa{x)xa{x) po{x) ~ / dxXa{x)po{x) giving the 
population of the metastable state a in the equilibrium distribution po{x). 

The time evolution of the probabilities Pa{t) now follows from the projection of the 
Fokker-Planck equation (1) onto the slow subspace, i.e. we assume the probability density 
to take the form (10) and determine the scalar product with Xa(x) on both sides of the 
Fokker-Planck equation. This yields 

a' 

Under the sum on the right hand side we expressed the coefficients Ca,i{t) by means of eq. (13) 
in terms of the probabilities Pa' (t) leading to the equilibrium weights Ua' in the denominator 
of the respective coefficient. For a ^ a' the coefficient (xa, Lxa'Po) /^a' coincides with the 



Rayleigh-quotient expressions for the rate from a' to a [17]. The integral is dominated 
by a neighborhood of the saddle point that connects the domains of attraction Va' and 
Va [19]. Within this neighborhood, for small noise, the localizing function Xa(x) can be 
approximated by the characteristic function of Da-, [20]. The resulting coefficient then 

coincides with the better known flux-over-population expression for the respective rate [3]. 
In particular, one can show that for sufficiently weak noise the expression [xa-i Lxa'Po) I'^o.' 
is positive for a ^ a' . Prom eq. (9) it follows that the sum of the coefficients over a vanishes. 
Hence, eq. (15) has the proper form of a master equation: 

d 



where 



— Pa(i) = ^ ra,a'Pa'{t) ~ ^ ra',a Pa{t) (16) 



iXa,LXa'PQ) 

ra,a' = (17) 



is the transition probability per unit time from q:' to a. 

The regime of validity of this approximation, of course, depends on the particular system 
and the degree of accuracy one would like to achieve. By its very nature this master equation 
is only appropriate for the long time dynamics of the Fokker-Planck equation and therefore 
requires a clear separation of time scales dividing the process into fast relaxations within 
the domains of attraction and slow transitions between the attractors. We note that various 
definitions of transition rates exist which are different from a physical point of view but 
which still yield the same result if a master equation provides a correct description of the 
long time dynamics. It turns out that the actual differences between these rate expressions 
are exponentially small in the noise strength, i.e. of the form exp {— Ay/ZcfiT} where AV 
is a relevant barrier hight and ksT the noise strength. So if we are ready to accept an error 
of, say 1%, all barriers in the system must be higher than approximately 4.5A;bT'. 



III. SLOW DRIVING 



After these preparatory considerations we turn to our main topic and consider a Fokker- 
Planck process of the kind of the previous section with the only difference that we now allow 
for slow changes of the parameters of the system. These could result from a slow drift of 
the environmental parameters of the system like e.g. the temperature, a slow increase of an 
external field, or a slow periodic force driving the system, to name but a few realizations. 
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We emphasize that the change of the parameters must be slow compared to the fast 
relaxation within each domain of attraction but no assumption is made how it compares 
to the slow intrinsic dynamics of the system describing the transitions between different 
domains of attraction. We assume that in the presence of these time dependent parameters 
one still can identify domains of attraction which do not merge or split as time changes but 
retain their identity. In other words, we restrict ourselves to changes of the parameters that 
do not lead to topological changes of the system. 

The resulting process is still Markovian and is governed by a Fokker-Planck equation: 

= L(t)p(t) (18) 

where the Fokker-Planck operator L[t) is of the form given in eq. (2) with time dependent 
drift K(x, and possibly also a time dependent diffusion matrix D(x, i). In particular, 
we assume that all parameter values which the system experiences in the course of time 
correspond to equilibrium systems, i.e. that for any frozen parameter value the system will 
reach an equilibrium distribution relative to which the system obeys detailed balance: 

L{t)po{t) = Pom+{t) (19) 

where the tilde denotes time reversal and po{t) is the multiplication operator with the equi- 
librium distribution po(x, t) that would be attained if the parameter values were frozen at 
the values that they assume at time t: 

L(Opo(x,0 = 0. (20) 

For vanishing noise the system will move toward the nearest attractor as defined by the 
momentary values of the parameters. During this relaxation time, the parameters of the 
system are supposed to change by such a small amount that they can be considered as kept 
constant. Only when the system has reached the attractor and stays there for a much longer 
time, the change of parameters will become noticeable and the system will follow the slow 
motion of the attractor as it results from the changing parameters. This means that the 
time rate of change of the deterministic drift is characterized by a small parameter e: 

|K(x,t)| 
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where is a typical relaxation time of the deterministic motion. Under this very condition 
many of the properties of the time homogeneous process carry over to the case with time 
dependent parameters. 

To each attractor a then there belongs a domain of attraction D^if) together with its 
characteristic function xi°''(x, t) which is 1 on T>a{t) and elsewhere in E. The domains of 
attractions also slowly depend on time and so do the characteristic functions. If the system 
is weakly perturbed by noise the characteristic functions will change into smooth functions 
Xai^, t) that coincide with the noiseless functions everywhere except in a thin boundary 
layer at &Da{t), where a steep but smooth transition from 1 to takes place. As in the time 
independent case, they arc solutions of the homogeneous backward equation that approach 
1 in the interior of T>a{t) and vanish on all other attractors a' ^ a. Further one also can 
represent the slow dynamics of the probability density in terms of linear combinations of 
Xa(x, t)po{K,t): 

m 

p(x, t)^Yl Ca(OXa(x, 0P0(x, t) (22) 

The probabihties to find the system in the state a is analogously defined as in the case with 
constant parameters: 

Pait) = {Xa{t),p{t)) (23) 

The time rate of change of Pa{t) consequently has two contributions, resulting from the 
derivative of Xa(x, t) and of po(x, t): 

The latter probabilistic contribution can be expressed in terms of the Fokker-Planck equation 
yielding 

(xa, ^p{t)] = J]r«,„'(i)Pa(i) (25) 

V 'a' 

where we used eq. (22) and introduced the time dependent transition rates ra,a'(^) in an 
analogous way as in the time homogeneous case, see eq. (17): 

r.,.Kt) = (26) 

Ua' \t) 

It determines the rate of change of the probability "Pak^^ caused by transitions from ol to a. 
Here n^' (i) is the population of the metastable state a' in the frozen equilibrium distribution 
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Po{x,t). It is given by: 



na{t) = {Xa{t),Xa{t)po{t)) ^ {Xa{t) , Po{t)) (27) 

where the term Xq,(x, t) is neglected in the second equation. This approximation holds up 
to exponentially small terms in the noise strength. 

The other contribution, {dxa{t)/dt,p{t)), describes the change of the probability caused 
by the geometric change of the domain T>a{t). It also is linear in the probabilities Pa'{t) and 
can be written as 

(^XaW,PW) = X]^a,a'(t)Pa'(^) (28) 

where 

The coefficients ga,a'{t) are proportional to the smallness parameter e because they contain 
the time derivative of Xai^i't)- The proportionality factor roughly is of the same order of 
magnitude as the transition rate ra,a'{t)'- 

^^^^ = e»(e) (30) 

The dependence of this ratio on the noise strength is discussed for a particular example in 
section V. For a sufficiently slow parameter change, i.e a small value of e, the geometric 
contribution to the time rate of change of the probability Pa{t) can be neglected and a master 
equation results that is completely determined by the instantaneous rates ra,a' (t) '■ 

^Pa(^) = ^oc,a'{t)Pa'{t) - ,a{t)Pa{t) (31) 

where we used that the time dependent rates ra,a'(t) (26) have the same general form and 
therefore the same formal properties as the time independent ones (17): Those for a ^ a' 
are positive and the sum over the first index vanishes: ra,a'{t) = 0. This master equation 
with time dependent rates is the central result of the present paper. For faster driving the 
geometric contributions of the rates have to be taken into account. However, they may 
become negative, and therefore formally negative rates may result if the driving is too fast. 
This indicates that then the instananeous eigenfunctions of the Fokker Planck equation no 
longer provide a good basis. 
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IV. DECORATING THE META-STABLE STATES 



Often the knowledge of the probabihties Pa{t) to find the system at time t in the meta- 
stable state a is not sufficient. For example one may be interested in the average position x 
of the system or its single- or multi-time statistical properties. We here show how the time 
dependence of these quantities on the slow time scales of the transitions and the external 
driving can be retrieved from exactly the same information that is necessary to determine 
the transition rates from the Fokker-Planck equation. For that purpose we first consider the 
average of an arbitrary function /(x) of the position: 

{fit)) = J r^f{^)p{^,t) 

= J2{f{t)\a)p^{t) (32) 

a 

where {f{t)\a) denotes the expectation value of f{x) under the condition that the system 
resides at t in the domain of attraction (t) . Using the long-time behavior of the probability 
density p{t) as given by eq. (22) one finds: 

{f{t)\a) = J dM(x)p(x,i|«) (33) 

where the conditional probability p{x,t\a) to find the system in the continuous state x e 
V^{t) is given by 

p(x,t|a) = ^^^Po(x,t) (34) 

Similarly, one finds for the correlation of two functions /(x) and g{'K) at times t and t + r 
being separated by a positive time r that is long compared to the fast time scale of relaxations 
within the domains of attraction the following expression: 

{fit + T)g{t)) = J2{fit + r)\a) {git)\a') p{a, t + T\a', t) p^.{t). (35) 

a,a' 

Here p{a,t-\-T\a\ t) is the conditional probability to find the system in the metastable state 
a at time t + t provided it was in the state a' at the earlier time time t. This conditional 
probability is the solution of the master equation (31) subject to the initial conditions 
p{a,t\a',t) = Sa,a'- Eq. (35) holds for arbitrary functions /(x) and g{x.)- Therefore the 
conditional probability of finding the system at time t + T at the continuous state x if it was 
at the earlier time t at y takes the form 

p(x, t + r|y, t) = ^ p(x, t + T\a) Xa'iy, t) p(a, t + r|a', t). (36) 

a,a' 
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For an independent derivation of tfiis result see the Appendix A. 

Yet another analogous way exists to characterize the continuous process on the long time 
scale: One takes the discrete process z(t) assuming the values a = 1 . . .m according to 
the master equation (31) and decorates the states with a random point in the state space, 
^{z{t),t), depending on time and the particular state a = z{t) that is realized at the time t. 
If the probability density of the decoration X(a, t) is chosen as the conditional probability 
density p{x.,t\a) given in eq. (34) the mean values and the correlation functions of the 
process X(x, i) coincide with the expressions (32) and (35), respectively, and accordingly 
are characterized by the conditional probability (36). The decorated process is Markovian 
and has the same single and two time probability density as the continuous process x(t) on 
the long time scale and hence coincides there with it: 

^(t)^X(z(t),t) (37) 



V. BISTABLE OVERDAMPED OSCILLATOR 



We here consider the example of an overdampcd bistable Brownian oscillator that is 

periodically driven at a frequency fl which is slow compared to the typical relaxation rates. 

Its dynamics is given by the Smoluchowski operator 

LU.9_9VM^^_,^ (38) 
dx dx dx"^ 

where is the noise strength which is proportional to the temperature 9 of the fluid 
surrounding the oscillator, and V{x,t) denotes the time-dependent potential: 

V{x, t) - ^x^ - ^x'^ - Ax sm{flt) . (39) 

Here A denotes the strength of the periodic force. Throughout dimensionless units are used. 
It is supposed to stay within the limits \A\ < 2/{3^/^. Actually it must keep some finite 
distance from these limits in order that the following asymptotic theory applies. Under this 
condition on the strength of the external force, the potential has three stationary points 
that are solutions of the algebraic equation: 

x^ - X = Asm{nt) (40) 

As time varies they form three branches: Xi{t) and X2{t) are those tracing the two local 
minima that at t — assume the values Xi{0) — —1 and X2(0) = 1, respectively, and Xb{t) 
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gives the location of the local potential maximum between the minima. At i = it is located 
at Xb{0) — 0. The corresponding extreme values of the potential V{x, t) are denoted by: 



Va{t) = V{xait),t) for a = 1,2 

Vb{t) = V{xh{t),t) (41) 
The instantaneous well and barrier frequencies are defined accordingly as: 
^a{t)^ yaMpZ^ V34(i) - 1 fora = l,2 

Mt) = 4^^^^ = ^Jl - (42) 

For later use we give the time derivatives of the barrier position and the barrier frequencies. 
For the position one finds from eq. (40): 



and for the frequency it follows: 



VtA 

Xb{t)^ ^cos(Qi) (43) 



Mt) = -3^x,(t) (44) 



The instantaneous stationary solution of the Smoluchowski equation Ls{t)po{x,t) = is 
given by the Boltzmann distribution: 



Po 



>,i) - Z-i(i)e-'^^(^'*) (45) 



where for small noise the partition function is given by 

^^^^^ [2^ f exp{-/3V,{t)} ^ exp{-/3^(t)} | ^^^^ 

Here algebraic correction terms of the order 0{l/{pAVa{t))) were neglected. The domains 
of attraction of the instantaneous locally stable states Xa{t), a = 1,2, extend from minus 
infinity to the instantaneous barrier and from there to infinity: T>i{t) = {— 00, Xb{t)), and 
V2{t) = {xb{t), 00). 

The corresponding localizing functions Xa{x, t) are the solutions of the backward equation 

ujb{tf{x - Xb{t))^ + P-'^ \ Xa{x, t) = (47) 
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with boundary conditions: 



X2{x,t) 



The solutions readily are found: 



1 for a; — >• — oo 

for x — > oo 

for X — > — oo 

1 for X — > oo 



(48) 



Xl{x,t) = hTfc{Ub{t)^/pj2{x - Xb{t))) 

X2{x,t) = l-xi{x,t) (49) 



In the present model there is only a single variable, the coordinate x which transforms 
evenly under time reversal, hence, Xai^, t) = Xa{x, t). Using the eq. (49), the scalar products 
iXait), Ls{t)xa'{'t)Po{'t)), 0(,0i' — 1,2 Can be expressed by a single one that we denote by 
q{t): 

{Xi{t),Ls{t)x2{t)po{t)) = {x2{t),Ls{t)xi{t)po{t))^ 
- ixiit), Ls{t)xi{t)po{t)) = - {X2{t), Ls{t)x2{t)po{t)) = q{t) (50) 



where q{t) takes the form: 



q{t) = J dx (^^^^) ' Po{x, t) (51) 



Using eq. (49) for weak noise it can further be simplified as 

= -mvw <^^' 

The populations na{t) of the metastable states a = 1, 2 in the frozen equilibrium distri- 
bution become at weak noise: 

Using eq. (26) this gives for the rates the instantaneous expressions [3]: 

nAt) - -r2At) - ^^^ai^e-Z^^^^W (54) 
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A. The geometric correction to the rate 

Now we come to the discussion of the corrections of the instantaneous rates that are 
determined by the geometric contributions {dxa{t)/9t,Xa'{t)Po{t)): o;, a' = 1,2. For the 
relative magnitude of the geometric correction to the instantaneous rate r2,i(t) one finds 
after some algebra: 



As expected, the error is proportional to the driving frequency Q via the time derivative 
Xb{t)- The remaining f2-dependence may be absorbed in a rescaled time r = fit. An analytic 
solution of the integral in eq. (55) is not known. In Fig. 2 the error divided by the driving 
frequency is shown as it results from a numerical evaluation for different values of the driving 
strength A and inverse temperature as a function of time. Both increasing driving strength 
A and inverse temperature /3 lead to an increase of the extrema of the error that must be 
compensated by a smaller driving frequency in order that the master equation with the 
instantaneous rates provides a valid description. 

As an average measure we introduce the root mean square error E — ij^ dte^{t)/T\ 
which is strictly proportional to the driving frequency Jl. The proportionality factor E/Q 
increases as a function of the inverse temperature first as a power law and changes to an 
exponential growth for large values of /3, see Fig 3. 

We propose the following physical explanation of this effect: At large f3 the noise is very 
small such that a particle that has reached at some time the top of the barrier has a very 
small chance to be pushed to either side by the noise. During the same time when it resides 
there, the potential moves such that the particle no longer will sit on top of the barrier. 
Depending on the direction of the motion of the potential the particle may now be on the 
opposite side of the barrier or on the original side where it came from. In the former case 
the geometric correction leads to an increase of the rate and in the latter case to a decrease 
of the transition rate. Because it then has again to overcome a barrier the increase of the 
rate is exponentially large in the inverse temperature. 



e{t) 



92,i{t) _ {dx2{t)/dt,Xiit)poit)) 
r2,i{t) {X2it),Ls{t)xi{t)po{t)) 




(55) 
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FIG. 2: The relative error of the rate oscihates as a function of time r = 0,t. It vanishes at 
the extrema of the driving force. Panel (a) shows e(r)/il for different inverse temperatures /3 = 
30, 60, 90 at the driving strength A = 0.15 and panel (b) for different driving strength A = 
0.1, 0.15, 0.2 at /3 = 60. The parameter values are indicated close to the respective curves. 




500 1000 1500 2000 

beta 



FIG. 3: The asymptotic exponential growth of the frequency independent ratio E/0, is shown for 
the driving strength A = 0.12. The straight line corresponds to a barrier height 0.97 x 10~'^. 
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Finally we note that the deviation from the frozen rate on the noise strength depends 
on the particular model. For example in symmetric systems with varying barrier height the 
deviation only grows with the square root of the inverse temperature rather than exponential 
as in the above case where both the location and the height of the barrier changes in time. 

B. Average motion 

We have projected the dynamics of a a process with multiple metastable states and 
slowly time dependent parameters onto the slow subspace defined by the eigenfunctions of 
Fokker-Planck operator 

As a simple example of the decoration of the two metastable states we come back to the 




-0.04 -0.02 0.02 0.04 
F(t) 

FIG. 4: The response of the average position on the driving force for the driving strength A = 0.05 
and driving frequency O = 10~^ strongly depends on the magnitude of the noise. For /? = 35 
the noise is relatively large such that the average position closely follows its adiabatic value and 
therefore shows almost no hysteresis; for /? = 48 it stays behind the force, and for /3 = 75 it is 
hardly influenced by the force. The hysteresis curves are traced in the mathematically positive 
sense. 
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bistable oscillator and consider the asymptotic motion of its average position. Prom eq. (32) 
it results as 

(xit)) = {x{t)\l)p,it) + {xit)\2)p,it) (56) 

where {x{t)\a) denotes the conditional average of the position in the well corresponding to 
the metastable state a: 

{x{t)\a) = / dx x ^"^^'!^ po{x,t) 

~ Xa{t) (57) 

In the last equation we disregarded small contributions of the order 0{l/{u;a{t)y/P)). The 
probabilities Pait) follow as the asymptotic solutions of the master equation: 

P,it) = \_,KiT) ^"^^*^+ / rf«e-(^W-^(^))r,,.(.) 

P2{t) = l-P2{t) (58) 

where 

K{t)= [ ds{n,,{t)+r2,i{t)) (59) 
Jo 

Rather than considering the time dependence we here show in Fig. 4 the dependence of 
the position on the driving force F{t) — Acos{Qt). In general the position lacks behind the 
force and therefore the dependence is hysteretic. Only if the driving is either very slow or 
very fast compared to the average rates there is no hysteretic behavior. In the slow case 
many transitions occur before a change of the potential becomes sensible. The asymptotic 
probabilities Pait) take the form of the stationary probabilities for the frozen rates. They 
then depend on the force but not on its rate of change. In the opposite limit the asymptotic 
probabilities see only the average rates and hence are time independent. 



VI. SUMMARY 



In this paper we studied the dynamics of externally driven systems with multiple 
metastable states in the semiadiabatic limit. The full dynamics of the system compris- 
ing the fast and the slow time scales of the intrawell relaxations and of the external driving 
as well as of the transitions between the metastable states, respectively, is supposed to be 
Markovian and continuous and, hence, described by a Fokker-Planck equation. We showed 
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that the long time dynamics of this Fokker-Planck equation is equivalent to the dynamics 
of the master equation for the transitions between the discrete metastable states. The rates 
that determine this master equation are time dependent taking values as if the parameters 
were frozen at their instantaneous values. This kind of kinetic description is not new and 
has long been used in the literature [13, 22]. The present work, however, provides three 
new aspects. First the master equation does not result from an educated guess. It rather 
is obtained as the result of the projection of the dynamics onto the slow subspace of the 
Fokker-Planck equation. This systematic approach additionally gives a quantitative crite- 
rion when deviations from the master equation with frozen rates must be expected. This 
criterion docs explicitly take into account the time change of the parameters and in contrast 
to other works [11] is not based on a comparison of time scales of the frozen dynamics with 
the rate of change of the parameters. As a particular result we found for the dynamics of 
a periodically driven overdamped Brownian oscillator that the maximal driving frequency 
depends on the noise strength and decreases with decreasing noise strength. Asymptotically 
this dependence even is exponential. Future work still has to show whether one can achieve 
an improvement of the dynamics if the so called geometric corrections to the rates are taken 
into account in cases when the time scale separation between the fast intrawell relaxation 
and the changes of the parameters becomes less pronounced. The third main outcome of 
the present approach is a partial retrieval of the underlying dynamics of the continuous vari- 
ables by means of a proper decoration of the discrete states with time dependent random 
variables. 

The necessary time scale separation forced us to exclude the occurrence of bifurcations 
of the deterministic dynamics caused by the change of the parameters. Such bifurcations 
are accompanied by a slowing down of the deterministic dynamics and at the same time by 
a lowering of a barrier height and consequently by a decreasing time scale of transitions. 
Hence, for time dependent parameters, the time scale separation between interwell and 
intrawell dynamics is violated within a time window around the instant of bifurcation. It 
seems plausible that it should be possible to cut out this time window and to bridge it 
by some connection condition for the probabilities of those states that are involved in the 
bifurcation. 

Apart from the separation of the time scales, another assumption was made in the present 
work. First, we assumed that for fixed values of the parameters the system reaches a state of 
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thermal equilibrium and hence it obeys detailed balance. We think that this assumption is 
not really essential for our results but it simplifies the analysis considerably. Without detailed 
balance more general limit sets of the dynamics can occur apart from limit points such as 
limit cycles and chaotic states. In general, the stationary densities are not known in absence 
of detailed balance, the weak noise asymptotics is plagued by notorious nonanalyticities 
[23], and rather little is known about transition rates in nonequilibrium systems. On the 
other hand there are many important systems that are driven out of equilibrium by time 
dependent parameters. 
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APPENDIX A: THE CONDITIONAL PROBABILITY AT LONG TIMES 

For a sufficiently long time lag r the conditional probability p{x,t + T\y,t) of the 
continuous process x{t) can be expressed in terms of the basis functions Xa{y,t) and 
Xa{x,t + T)po{x,t + r) spanning the slow subspaces of the backward operator at time t 
and the forward operator at time t + t, respectively. One therefore can write: 

p{x, t + T\y, t) = J2 P'; t, r)p{x, t + r|/?)x/3'(y, t) (Al) 

where we expressed the basis function Xa{x,t + T)po{x,t + r) by p{x,t + T\a)na{t + r), 
see eq. (34) and introduced the yet undetermined coefficients d{a, a'; t, r) into which the 
time dependent factor na{t) is absorbed. We now multiply the conditional probability of 
the continuous process p{x, t + T\y,t) with p{y, t\a'), integrate over the domain of attraction 
1^a{t+T) and obtain the conditional probability of the discrete process p{a, t+T\a', t). Using 
the ansatz (Al) for p{x, t+r|t) we then find for the coefficients d{a, a'; t, r) = p{a, t+r|, a', t). 
Together with eq. (Al) this yields the expression eq. (36) for the conditional probability of 
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the continuous process as claimed in Sect. IV 
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